Crumpled Wires and Liouville Field Theory 
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The study of the packing of a length of wire in a two dimensional domain is done using techniques 
of conformal maps. The resulting scaling properties are derived through the Coulomb gas formalism 
of Conformal Field Theories. An analogy is drawn between this system and the 2D gravitational - 
or Liouville - coupling of non-conformally invariant matter in two dimensions. 
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Introduction — Crumpling is a general problem in 
Physics in a wide range of phenomena including mate- 
rials science, metal foils and organic matter [1], which 
poses interesting challenges in both theoretical and ap- 
plied perspectives. Recently some attention has been 
drawn to the problem of crumpling of wires, especially 
in two dimensions. Relevant applications include pack- 
ing of DNA strands in viruses, general polymer chains 
and confinement of elastic rods [2] . 

Ideas of scale invariancc have been pervasive in criti- 
cal phenomena since its inception, and in two dimensions 
its natural conclusion, conformal field theory, has been a 
valuable tool since [3], especially in the context of ran- 
dom surfaces and integrable models. These models also 
provide a concrete realization of two dimensional quan- 
tum gravity, in its representation as a Liouville field the- 
ory [4] . Integrability via Toda hierarchy has exposed an 
exact equivalence of 2D quantum gravity and the prob- 
lem of Laplacian growth [5] . Given the lack of intuitive 
understanding of quantum gravity, it would be desirable 
to find more amenable problems which possess analogies 
with more realistic gravity problems, such as the coupling 
to non-conformally invariant matter. 

In this letter we show an equivalence between the prob- 
lem of two-dimensional packing of wires and Liouville 
field theory coupled to matter in an "almost conformal" 
manner. The breaking of conformal symmetry stems 
from the boundary conditions and are hence tractable 
by conformal perturbation theory. We also give a simple 
scaling argument that predicts a phase transition where 
the length of packed wire is a maximum and the critical 
exponent of the energy. 

The set up — Let D be the interior of the container, 
seen as a region of the complex plane. For simplicity, 
take D to be the unit disk. The packing process begins 
with a diameter of the disk being continuously deformed 
so as to increase its length while keeping the angles at 
the extremities fixed, which is called an injection. At a 
given wire length L, consider z^it), the map that takes 
the initial diameter t G [—1, 1] to the configuration of the 
injected wire C. We will want to view the function 2l(0 
as the restriction of a transformation (or diffeomorphism) 
of the whole domain D into itself. We will argue below 



that the condition of non-self-intersections of the curve C 
can be cast in terms of the diffemorphism of the domain 
D by saying that the regions above and below the curve 
C, -D+ and _D_, are respectively mapped to the upper 
and lower semi-disc by analytic functions. 

Such extension of ZL{t) to the whole domain D pre- 
serves its manifold structure and as such the function 
ZL{t) can be seen as locally analytic. In the next section 
we will further explore this property. Also, the restric- 
tion of such maps to the boundary of D, the unit circle 
in our case, forms a group, DiffS^ with ubiquitous 
applications to conformal field theory. Seeing D as a 
complex manifold, the restriction enforces that we can 
cover the image of D with an atlas consisting of two sim- 
ple analytical maps, with an ovclap region that contains 
a neighborhood of ZL{t). By definition, a simple map 
takes simple curves to simple curves and therefore the 
restriction forbids intersections [6]. The transition func- 
tions defined in the overlap will themselves be analytic in 
order to warrant non-self intersection. The construction 
will be explicited in the last section. 

Physically, the importance of this restriction arises as 
one injects wire into the domain, or as L increases. As it 
happens, the system tries to minimize the energy creating 
points of self-intersection. The transition is actually akin 
to the generation of vortices in the XY model, in which 
the angle of the tangent to the wire 6{t) goes over a full 
loop before reaching the other side of the domain 9{t = 
1) =e{t = -I) + 211. 

Let us now consider the map z{w), with L undeter- 
mined. The curvature k of the curve zit) and its element 
of length M are given by 
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where z{t) = z*{t) is the complex conjugate of the map. 
The problem of adiabatic growth is thus stated: let z(w) 
be the restriction of the diffeomorphism of the unit disk 
D to D defined in a neighborhood of w real. The con- 
formation of the wire inside the domain z(t) is such that 
the free action 
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is a minimum. The Lagrange multiplier fi implements 
the constraint that the minimization has to be done while 
keeping the total length fixed: 
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L. 
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At any new injection of wire into the system, it is then 
assumed that there is enough time for the elastic energy, 
given by the first term of F, to be minimized. In prac- 
tice, this time can be of the order of several seconds for 
ordinary materials. The general form (2) is invariant by 
reparametrizations of the curve t t'{t) and as such it 
will be reduced to the usual forms by a suitable choice 
of the parameter. Moreover, it will yield, via a suit- 
able change of field variables, to the relation between the 
problem of crumpling and Liouville field theory. 

The free energy (2) is easily recognizable in the proper 
length parametrization z = e'^ , with 9 being the angle of 
the curve with a fixed direction: 
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One can then fix the final position of the wire by means 
of a Lagrange multiplier, arriving at an expression for the 
free energy: 
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Since the constraint is global, one can treat the Lagrange 
multiplier as a constant coupling. This form is unsuitable 
to treat the various other constraints, but it is an ade- 
quate approximation when the free energy is dominated 
by the elastic bending. The value of A is determined in 
terms of L by imposing that the endpoints are separated 
by a diameter At = 2. 

The Schwarz function and Liouville Theory — 
Another interesting parametrization of (2) la t — z, valid 
on the curve C . In this case, the relevant map is z = 
S'(z), where S{z) is called the Schwarz Function [7]. The 
tangent vector to the curve z(t) is given by the relation: 



dzdz = S'{z)dz'^ 
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whereas the curvature is given by 

de _^ i dS' _ i dz dS' _ i S" 
di ~ 2S'~d£ ~ 2S'de~d^ ~ 2 {S')^/^ 
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One can then try to minimize the functional (2), but it 
is easier to consider (5), where the constraint of constant 
length is already implemented. From it results the re- 
lation ^ + i^ie:'^ - e^''^) = 0, which in terms of the 
Schwarz function is 
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where {S; z} = S"'{z)/S'{z) - i/2{S"{z)/ S'{z)f = T is 
the Schwarzian derivative of S'(z). Upon the substitution 
S"(z) = e'^'''^^^\ the equation can be written as 



T + V 



=2701 



(9) 



which is recognized as the (chiral) "gravitational Ward 
identity" of Liouville Field theory [8] coupled to "matter 
field" whose stress energy tensor is given by the expo- 
nential terms. The constant 7 will be fixed using scaling 
arguments below. 

Strictly speaking, (9) is valid only on the curve C, 
but the requisite of non-self intersection of the curve 
means that the equation above should really be thought 
of as a analytical continuation of the function (j){t) to a 
neighborhood of C. By the discussion of the preceed- 
ing section, this in turn allows us to continue the func- 
tion, and S{z), into either chart of the domain D. In 
fact, the Schwarz function does allow for a decomposi- 
tion S{z) = S^{z) + S-{z) such that S±{z) is analytic 
at either side of C [7]. One can further enforce that 
S±{z) thus constructed keep the circle invariant, but for 
now it suffices to say that the condition of simplicity is 
equivalent to the non-vanishing of S'j. and hence that (jj 
is well-defined. 

We will be interested in the scaling properties of ob- 
servables like the elastic energy (4): 

\2 
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and the total length 
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(11) 



dz^/S' = 
ic JC 

We note that since C belongs to both charts, one can 
then deform the contour to "half a perimeter" of S^, any 
line connecting w ^ ~1 to w ~ 1, without changing the 
value of either observable. 

Each of those observables has a definite conformal di- 
mension, which in this case is offset by the "improve- 
ment term" ^d^cfi in (9). This is best dealt in the so- 
called Couloml) gas formalism for Conformal Field The- 
ories (CFTs) [9]. Given the similarity between T in (9) 
and the free boson, let us start without the improvement 
term or the potential. For a free boson the stress energy 
is To = — and the Operator Product Expansion 

(OPE) of correlation functions are computed from the 
Green's function {(j>{z)(f>{z')) — — log(z — z'). One can 
then compute the OPEs between T and the type of ob- 



servables relevant to us Vp{z) 
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2{z-z'y . . 

where "reg." stands for a regular part which implement 
the effect of boundary conditions. Since To (2) is the gen- 
erator of conformal transformations, the constants in the 
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leading order term gives us the behavior of Vp{z') under 
scale transformations, being in fact the operator's scaling 
dimensions. We then have A[V^] = — /3^7^/2. Another 
type of observable important for the discussion is a de- 
scendant, which comes about by the application of the 
Fourier components of Tq{z) to Vp: 

L„v^.H-2^£(^^n^)^.H. (13) 

Because of the Virasoro algebra, such fields have dimen- 
sion A[yg] + 1. We will be particularly interested in 
LiVp = f3"f{d(l))eP^'>' . 

The usual interpretation [4] of the improvement term 
in (9) is that the otherwise free conformal field is really 
coupled to the curvature of the base space (the w plane 
in our case), via the term (p^/gR. The system is no longer 
invariant over fixed translations of 0, but the variation 
is proportional to the Gauss-Bonnet term J yJgR = Attx 
and hence a constant for fixed topology. This extra term 
will contribute to the stress energy tensor since the lat- 
ter is just the variation of the action with respect to 
the metric, and y/gR does have a non-trivial variation 
SiVaR) = M9'''9'"'-9''''9''')'^aVb59cd- Integrating by 
parts yields a contribution by a term proportional to the 
second derivative of 0, Qd'^cj), just like the one appearing 
in (9), where Q = l/2j. 

Owing to the extra term in (9), the OPEs between 
T(z) and the local operators Vg = e'^^'^ will be shifted, 
and their new conformal dimensions will be 

A[Vp] = QP^-^ = ^{l-P^'). (14) 

As far as the descendants arc concerned the effect of the 
improvement term is in general non-trivial, but a very 
fortunate coincidence emerges for the case we are inter- 
ested in: the first Fourier coefficient of TVp is not affected 
by its presence. One can verify this by expanding the "ki- 
netic" , or "elastic" part of the stress energy in terms of 
Fourier components of the field (j>. In short, the rela- 
tion between the scaling dimensions of the field and its 
first descendant remains valid even when the improve- 
ment term is introduced. 

We are now in position to fix 7 by requiring that both 
potential terms in (9) have the same scaling dimensions, 
and substituting values one has 

1-72 3 - 97^ 

A[^i] = ^^, A[V^]^^^, (15) 

equating these gives us 7 = ±1/2 and A[Vi] = A[V3] = 
3/8. The scaling properties of the relevant observables 
E (10) and L (11) are dealt with by different means. L 
is the integral of a operator of dimension A[Vi] = 3/8, 
but since it is integrated over a path, this number will 
be reduced by one as far as scaling transformations are 
concerned: A[L] = —5/8. Now the integrand in (10) is 



the square of the first descedant of the operator Its 
dimension is then given by A[E] = 2(A[y_i/2] -I- 1) — 1 = 
7/16, which means that these two operators have well 
defined scaling relations with respect to a change of the 
parameters, like A in the Ward identity (9). As there 
is only one dimensionful parameter on the theory, the 
constant a, one can summarize both behaviors in the 
single equation, stating how the energy scales with L: 

E^\L,,it-L\-i>, (16) 

where Lcrit is a function of a and the particular form 
of the domain D. This hints at a phase transition in 
which there is a maximum length of wire that can be 
injected in the domain. One also notes that given the 
scaling properties of icrit, it should increase with the 
size of the domain, and it follows that the energy will 
become more and more concentrated as the size of the 
domain increases. This is a feature the system shares 
with gravity. One should note that the exponent above 
should be valid as long as the system is in the elastic 
regime. 

The Boundary conditions — The considerations 
made thus far are valid regardless of the boundary con- 
ditions. Although the details of the constraint set the 
particular form of the potential in (9), this is not partic- 
ularly relevant for the scaling properties studied in the 
last section. Here we will present a procedure of how to 
deal with it from the analytical map standpoint. 

Let us start from the consideration of the map z(u'). 
Since it takes the unit disk to itself, the map ip{'w) = 
— ilog2;('u;) is real on the unit circle. It is related to 
the function (l>{w) = ~i\ogdS{w) by (j){w) = ip{w) — 
i\ogdip{w) + ^. One would then like to impose on 0(w) 
constraints that reflect 'ip{'w = 1/w) = tp{w). 

The constraints are conveniently written in terms of 
the harmonic functions u = ((/>(w) -I- 0(w))/2 and v = 
{(j){w) — (j){'w))/2i, conjugate to each other. On the unit 
circle, we then have 

u(w,w ^ l/w) ^ 'ijj{w) + ilogw ^ 
viw^w = 1/w) = —log{iwdTp{w)). 

The condition that ip{w), and then z{w), is locally ana- 
lytic in the disk is that dtp{w) ^ on either chart. This is 
tantamount to saying that viw^w = 1/w) docs not have 
singularities on cither chart. Combining both equations 
above one has e^'" = —dtu + 1, where dt is the tangen- 
tial derivative to the boundary. In turn, d^u is equal to 
minus the normal derivative of v, since they are harmon- 
ically conjugate. It is also clear that the numerical value 
1 represents the curvature of the boundary. Hence the 
constraints satisfied by the function v[w, w) can then be 
accounted for by the action: 

y=^( d^zdvBv + ^l di{vK + €-")■ (18) 
27r Jd 47r Jqjj 
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where K is the curvature of the boundary, constant in our 
case. The area element is (fiz = ^dz Adz. The Neumann 
boundary condition that ensues from the variation of (18) 
is exactly the one that arises in the correspondence with 
the sum over triangulated surfaces, a discrete version of 
two-dimensional quantum gravity [10]. 

The minimization of the elastic energy at constant wire 
length is given by the addition of the free energy (2). 
Writing it in terms of u and v, one has 

/I pi 
dt{ufe"+n j dte-^\ (19) 

We note that ^ is exactly of the form of a point particle 
moving in a + 1-dimensional space-time. As such the 
full action + ^ has the interpretation of a point par- 
ticle of mass [I interacting with two dimensional gravity. 
This action is also a natural candidate for branc dynam- 
ics in non-critical string theory [11]. For our purposes, 
the condition -^u = -£^v can be seen as a gauge choice of 
parametrization of the worldline of the particle. Unique- 
ness of this parametrization implies that the worldline of 
the particle is a simple curve inside the domain. We note 
that the addition of ^ means that the function v{w, iJu) 
defined above is no longer harmonic, meaning the break- 
age of the single analytical chart parametrization of the 
whole disk. Since the function is still harmonic at D±, 
one can define analytic maps in these domains which thus 
will comprise the atlas of D as alluded in the second Sec- 
tion. 

From the studies of Liouville field theory on a disk, 
the action y + ^ can be alternately interpreted as a 
deformation of pure Liouville by means of vertex opera- 
tors defined by (19). These are exactly E and L defined 
in (10) and (11), albeit in a different parametrization. 
It also follows that they satisfy the same OPE with the 
stress-energy tensor as computed above and hence the 
critical exponent should hold. The effect of the bound- 
ary is just to implement a identification between vertex 
operators Vp and Vg-p [12]. The two terms that com- 
prise the potential in (9), implementing the constraint 
that z{w = ±1) = ±1, are exactly of this form. The 
effect of the boundary is then to render one of them re- 
dundant, and thus the last ingredient of the action is 
the addition of the term A / d^ze""'"™. This constraint 
can nevertheless be dealt with by considering correlation 
functions of fields inserted at w = ±1. 

Given all the similarities explained above, one can thus 
posit that the transition for large wire length L — > Lq can 
be understood in terms of the OPEs of vertex operators 
(see [13] for a survey). In fact, since in this regime the do- 
main is covered in teardrop-like cells, one can conjecture 



that the system will undergo a Kosterlitz-Thouless type 
of transition by coalescence of "semi- vortices" (or vortices 
dipoles) where 9{z — ie) ^ 9{z + ie) + tt, corresponding to 
the point where the cell touches itself. Another interest- 
ing prospect is the interpretation of the action as gravity. 
It is believed that the non-conformal coupling may help 
gravity fluctuations that are important for inflation and 
dark energy [14], thus we hope to address this point in 
the simplified model proposed here in future work. 
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